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Nonthermal Emission from Star-Forming 
Galaxies 

Yoel Rephaeli and Massimo Persic 



Abstract The detections of high-energy 7-ray emission from the nearby starburst 
galaxies M82 & NGC253, and other local group galaxies, broaden our knowl- 
edge of star-driven nonthermal processes and phenomena in non-AGN star-forming 
galaxies. We review basic aspects of the related processes and their modeling in 
r~| ■ starburst galaxies. Since these processes involve both energetic electrons and pro- 

Qh! tons accelerated by SN shocks, their respective radiative yields can be used to ex- 

plore the SN-particle-radiation connection. Specifically, the relation between SN 
^ [ activity, energetic particles, and their radiative yields, is assessed through respec- 

c/5 ■ tive measures of the particle energy density in several star-forming galaxies. The 

^, I deduced energy densities range from i^( 10^ ' ) e V cm^^ in very quiet environments 

to 1^(10^) eV cm^^ in regions with very high star-formation rates. 



1 Introduction 



High star formation (SF) and supernova (SN) rates in starburst (SB) galaxies (SBGs) 
(^ [ boost the density of energetic nonthermal particles, whose main constituents are 

^f) • protons and electrons. Coulomb, synchrotron and Compton energy losses by the 

electrons, and the decay of pions following their production in energetic proton 

interactions with protons in the ambient gas, result in emission over the full electro- 

. magnetic spectrum, from radio to very high-energy (VHE, >100GeV) y-rays. The 

X' 

."?.■ Yoel Rephaeli 

School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 

and 

Center for Astrophysics and Space Sciences, University of California, San Diego, La JoUa, CA 

92093-0424, USA, e-mail: |yoelr@wise . tau . ac . il) 

Massimo Persic 

INAF/Osservatorio Astronomico di Trieste and INFN-Trieste, via G.B.Tiepolo 11, 1-34143 Trieste, 

Italy e-mail: |persic@oats . inaf . it| 



2 Yoel Rephaeli and Massimo Persic 

relatively high intensity emission in SBGs, as compared with emission form 'nor- 
mal' star-forming galaxies (SFGs), makes nearby members of this class the most 
likely non-AGN targests for 7-ray telescopes, such as Fermi and the (Cherenkov 
arrays) H.E.S.S., MAGIC, and VERITAS. 

Interest in 7-ray emission from SFGs clearly stems from the prospects for im- 
proved understanding of the origin and propagation mode of energetic electrons and 
protons and their coupling to interstellar media. This interest has been enhanced by 
recent detections of the two nearby SBGs M82 & NGC253 by Fermi (Abdo et al. 
2010a) and, respectively, by H.E.S.S ((Acciari et al. 2009) and VERITAS (Acero et 
al. 2009). M31, the closest normal spiral galaxy, was also detected by Fermi (Abdo 
etal. 2010b). 

A realistic estimate of the expected 7-ray emission requires a detailed account of 
all relevant energy loss processes of energetic electrons and protons as they move 
out from the central SB source region into the outer galactic disk. Calculations of 
the predicted X- 7-ray spectra of nearby galaxies were made long ago with varying 
degree of detail (e.g., Goldshmidt & Rephaeli 1995, Paglione et al. 1996, Romero 
& Torres 2003, Domingo-Santamaria & Torres 2005). A more quantitative numer- 
ical approach was initiated by Arieli & Rephaeli (2007, unpublished), who used a 
modified version of the GALPROP code (Moskalenko & Strong 1998, Moskalenko 
et al. 2003) to solve the Fokker-Planck diffusion-convection equation (e.g., Lerche 
& Schlickeiser 1982) in 3D with given source distribution and boundary condi- 
tions for electrons and protons. This numerical treatment was implemented to pre- 
dict the high-energy spectra of the two nearby galaxies M82 (Persic, Rephaeli, & 
Arieli 2008, hereafter PRA) and NGC253 (Rephaeli, Arieli, & Persic 2010, here- 
after RAP). The predictions made in these papers agree well with observations made 
with Fermi and TeV arrays, as will be discussed in the next section. 

Particle acceleration and propagation in galactic environments are largely similar 
in all SFGs. What mainly distinguishes a SBG from a normal SFG is the dominance 
of a relatively small central region of intense star formation activity. The overall 
validity of the numerical treatments of the two nearby SBGs provides a solid basis 
for generalizing the model to SFGs in general. 

We briefly review the calculation of steady-state particle spectra and their pre- 
dicted radiative spectra for the above two nearby SBGs, and discuss similar calcu- 
lations for conditions in a SFG. The particle energy density can be determined in 
several different ways. In order to assess and gauge the SN-energetic particle con- 
nection we compare estimates of the energetic proton (which dominate the) energy 
density in SFGs by three different methods, finding overall agreement, which pro- 
vides further evidence for the validity of the basic approach. 



2 Particle and Radiation Spetra in Starburst Galaxies 

Acceleration in SN shocks by the first-order Fermi process yields a power-law dis- 
tribution with index q>2 (e.g., Protheroe & Clay 2004) in a very wide energy range. 
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from a value close to the mean thermal energy of the gas particles (in non-relativistic 
shocks) to a very high value (>10''* eV). The accelerated proton-to-electron (p/e) 
density ratio, Np/Ng, in the source (either the SB or the full disk) region can be calcu- 
lated assuming charge neutrality (Bell 1978, Schlickeiser 2002). This ratio reaches 
its maximum value, (m^/msyi^^'''^ (for q>l; nie and nip are the electron and pro- 
ton masses), over most of the relevant range of particle energies, £>! GeV. (For the 
dependence of this ratio on particle energy, and more discussion on this and other 
relevant physical processes, see PRA and references therein.) 

The electron density in the source region, A^e, is inferred from radio measure- 
ments (of the same region); by adopting the theoretically expected expression 
Np/Ne = {mp/mgY'>^^''^, the proton density A^p can be deduced. The fit to the radio 
data provides both the normalization of the electron spectrum and the actual value of 
q, which is found to be somewhat larger than 2, even in the central SB region. In this 
procedure the electron population is composed of both primary and secondary elec- 
trons, with the latter self-consistently determined by accounting for the pion yield 
of energetic protons with protons in the gas. We note that the theoretically predicted 
value of the density ratio is valid in the source region, where energy equipartition 
is more likely to be attained since the relevant processes couple particles and fields 
more effectively than in the rest of galactic disk. 

The particle spectral distributions evolve differently as they propagate out from 
their acceleration region. Typically, the electron spectrum is most directly deduced 
from measurements of synchrotron radio emission. The inferred spectrum can be 
related to the source spectrum through a solution of the kinetic equation describing 
the propagation modes and energy losses by electrons and protons as they move 
out from their acceleration region. A very useful detailed description of the time- 
dependent spectro-spatial distribtion of protons, diffusing out of a region with a 
discrete population of acceleration sites, was recently given by Torres et al. (2012). 
This study elucidates the explicit dependence of the distribution on distance from 
the acceleration site, energy loss time, and the diffusion coefficient. It also follows 
the temporal evolution of the distribution towards a steady state. 

Since the estimated duration of a SB phase is ~ 10** yr, a timescale which is much 
longer than any of the relevant energy loss or propagation timescales for electrons 
and protons, in (essentially) all treatments a steady state is assumed to be attained. 
Since the calculation of particle steady-state spectra requires inclusion of all the 
important energy loss mechanisms and modes of propagation, the treatment is nec- 
essarily numerical. We have employed the code of Arieli & Rephaeli (2007), which 
is based on a modified version of the GALPROP code (Moskalenko & Strong 1998, 
Moskalenko et al. 2003), to solve the kinetic equation for Ni{Y,R,z), where / = e,p, 
7 is the Lorentz factor, and R and z are the 2D spatial radius and the coordinate 
perpendicular to the galactic plane. The exact Fokker-Planck diffusion-convection 
equation (e.g., Lerche & Schlickeiser 1982) was solved in 3D with given source dis- 
tribution and boundary conditions for electrons and protons. In addition to diffusion 
with an energy dependent coefficient, particles are assumed to be convected by a 
galactic wind with spatially varying velocity. 
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The dominant energy losses of high-energy electrons are synchrotron emission 
and Compton scattering by the FIR and optical radiation fields; these processes 
(obviously) depend on the mean strength of the magnetic field, B, and the energy 
density of the radiation fields, respectively. At energies below few hundred MeV, 
electrons lose energy mostly by Coulomb interactions with gas particles. At low 
energies proton losses are dominated by Coulomb interactions with gas particles. 
Protons with kinetic energy above the (range of) pion masses (~^140 MeV) lose 
energy mainly through interactions with ambient protons, yielding neutral (tt**) and 
charged (tt*) pions. Neutral pions decay into photons, while decays of tt* result in 
energetic e* and neutrinos. The proton and (total) electron components are coupled 
through the production of secondary electrons in n^ decay (following their creation 
in pp interactions). 

Measured synchrotron radio spectra provide the critically important information 
on the particle spectra and their overall normalization: Fitting the predicted radio 
emission to measurements fixes normalization of the steady state electron and - 
based on a theoretical prediction - proton energy distributions. From these mea- 
surements alone the electron density and mean magnetic field cannot be separately 
determined. To do so it is usually assumed that particle and magnetic field energy 
densities are equipartitioned. In our numerical treatment this approach necessitates 
an iterative procedure to solve for A'e, A^p, and the field strength at the center, Bq, 
given a measured value of the radio flux. 

Particles diffuse and are convected out of their source region. Diffusion is likely 
to be random walk against magnetic field inhomogeneities, with an estimated value 
of ~ 3 X 10^** cm"/s for the central diffusion coefficient. Convection is by galac- 
itc wind with a typical velocity of ^ 500 km s^^ (Strickland et al. 1997) in the 
source region. Based on Galactic cosmic-ray MHD wind models, we assume that 
the convection velocity increases linearly with distance from the disk plane (e.g. 
Zirakashvili et al. 1996). 

The other quantities needed to calculate the steady state distributions of electrons 
and protons are the densities of neutral and ionized gas in the central SB region and 
throughout the disk, the central value of the (mean) magnetic field and its spatial 
profile across the disk, and energy densities of ambient radiation fields (including 
the CMB). As discussed in PRA, it is assumed that magnetic flux is conserved in 
the IS ionized gas, so that the mean strength of the field can be related to the local 

9/3 

ionized gas density, «e, using the scaling B °^ rQ (Rephaeli 1988). If instead energy 
equipartition is assumed, and the magnetic energy density is scaled to the thermal 

1/9 

gas energy density, then the proportionality relation i?, B ^ n^.' " .In our work we have 
taken the ionized gas density profile to he. rig °^ exp{—z/zQ)/ ( 1 + {R/RoY), typically 
with Rq = \.5 kpc, and zq = 0.5 kpc (as deduced for NGC 253 by Strickland et al. 
2002). 

Given the measured radio fluxes from the central and full disk regions of the two 
nearby SBGs M82 and NGC253, shown in Figure 1 for the latter galaxy, and values 
of all the above quantities, the steady state particle spectra and their radiative yields 
were calculated using the modified GALPROP code. Here we present the results of 



Nonthermal Emission from Star-Fonning Galaxies 




0.5 

log(GHz) 



Fig. 1 Spectral fits to radio measurements of the SB and entire disk regions of NGC 253 (RAP). 
The solid line is a fit to the emission from the SB region; the dashed line is a fit to the emission 
from the entire disk. Data are from Klein et al. (1983, black dots), Carilli (1996, blue squares), and 
Heesen et al. (2008, green circles). 



this work; more details on the method and values of the input parameters can be 
found in PRA and RAP. 
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Fig. 2 Primary proton (dashed line), primary electron (solid line), and secondary electron (dashed- 
dotted line) spectral steady state density distributions in the central SB region NGC 253 (RAP). 



The high-energy photon spectra of NGC253 are shown in Figure 3 (from RAP). 
Emission levels depend mostly on values of the proton to electron ratio in the source 
region, on the magnetic field, gas density, and their spatial profiles. The basic nor- 
malization of the electron density is provided by the measured radio emission in 
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the source region; the variation of the electron spectrum across the disk is largely 
determined by synchrotron losses. Uncertainty in the estimated level of emission is 
largely due to the steep dependence of the electron density on the field. As argued 
by RAP, the central value of the (mean) magnetic field is unlikely to be appreciably 
higher than the value deduced in N253 (and also in M82), Bq ^ 200 ;UG. A lower 
field value would result in a reduced proton density and a lower rate of n^ decays. 
For a given radio flux the electron density would have to be correspondingly higher, 
resulting in higher bremsstrahlung and Compton yields, roughly compensating for 
the lower level of hadronic emission. Emission from n'^ decay depends linearly on 
the ambient proton density in the central disk region. 




Fig. 3 High-energy emission from the disk of NGC 253. Radiative yields are from Compton scat- 
tering off the FIR radiation field (dotted line), electron bremsstrahlung off ambient protons (dashed 
line), ;r'' decay (dashed-dotted line), and their sum (solid line). 



Detailed modeling of high-energy emission from M82 by PRA preceded its de- 
tection by VERITAS (Acciari et al. 2009) and Fermi (Abdo et al. 2010a). These 
observations are shown in Figure 4 (from Abdo et al. 2010a) , together with theo- 
retical predictions including that of PRA. The detected flux level agrees well with 
that predicted by PRA and by de Cea et al. 2009. Results from our similar treat- 
ment of the steady-state particle and radiation spectra of NGC253 were presented in 
RAP. The estimated high energy fluxes for this galaxy (RPA, Paglione et al. 1996, 
Domingo-Santamaria & Torres 2005) were also in the range measured by Fermi 
(Abdo et al. 2010a), and - in the TeV region - by the HESS telescope (Acero et al. 
2009), given the substantial observational and modeling uncertainties. These obser- 
vational results and predicted spectra are shown in Figure 5 (adopted from Abdo et 
al. 2010a); our integrated spectrum and estimated 1 a uncertainty region are marked 
by the dashed blue lines (inserted into the original figure). These two SBGs are the 
only extragalactic non-AGN sources that were detected at both GeV and TeV ener- 
gies. The FIR luminosities of these galaxies are at or somewhat below the nominal 
level for SBGs; clearly, their apparent high brightness is due to their proximity. 
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Fig. 4 High-energy emission from M82. VERITAS (Acciari et al. 2009) and Fermi measurements 
are shown together with predicted spectra. Figure is reproduced from Abdo et al. (2010a). 
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Fig. 5 High-energy emission from the disk of NGC253. Fermi and H.E.S.S measurements are 
shown together with predicted spectra. The figure is reproduced from Abdo et al. (2010b); the 
added dashed blue lines show our predicted spectrum with the estimated range of 1 CT uncertainty. 



A few other local and nearby galaxies were also detected at high energies: LMC 
(Abdo et al. 2010c), SMC (Abdo et al. 2010d), Andromeda (M31; Abdo et al. 
2010b), and the composite Sy2/SB galaxies NGC 1068 and NGC 4945 (Lenain et al. 
2010). The emission in these galaxies is mostly hadronic 7-ray emission, except for 
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NGC 1068 where emission from the active nucleus may be dominant. Among galax- 
ies whose detected high energy emission is powered by stellar activity, Arp 220 has 
the highest SFR in the local universe, yet an attempt to detect it by MAGIC (Albert 
et al. 2007) was not successful due to its relatively large distance. 



3 Estimates of Cosmic-Ray Energy Density 

Active star formation in galaxies leads to acceleration of protons and electrons via 
the Fermi-1 diffusive shock acceleration mechanism in SN remnants. Under equilib- 
rium conditions in a galaxy, a minimum-energy configuration of the magnetic field 
and the energetic particles may be attained. Energy densities of particles and mag- 
netic fields may then be in approximate equipartition, implying that the energetic 
proton energy density. Up, can be deduced from the detected level of synchrotron 
radio emission. In this radio-based approach Up can be estimated if the source size, 
distance, radio flux, and radio spectral index are known. 

In a y-based approach. Up can be obtained from the measured GeV-TeV spectral 
flux, which is mostly due to p-p interactions, as described the previous section. Only 
recently have such measurements become possible, at present only for 10 sources. 

In the SN method, with an assumed fraction of SN kinetic energy that is chan- 
neled into particle acceleration. Up can be estimated if the size of the star-forming re- 
gion and SN rate are known, as well as an estimate for the proton residence timescale 
from the presence (or absence) of a galactic wind emanating from the star-forming 
region. This timescale is largely determined by the advection timescale {^ lO^yr) in 
SBGs, and by the Tr'^-decay timescale ('-^ 10^ yr) in low-SFR (quiescent) galaxies. 

Expanding on our previous work (Persic & Rephaeli 2010), we show that the 
three methods give consistent results for Up for a sample of 10 galaxies with widely 
varying levels of star formation activity, from very quiescent to extreme SBs. These 
are the only galaxies of their kind for which y-ray data, in addition to radio data and 
SN rates, are available (see Table 1). 



3.1 Particles and magnetic field 

The population of NT electrons consists of primary (directly accelerated) and sec- 
ondary (produced via n"^ decays) electrons. While the exact form of the steady-state 
electron energy spectrum is not a single power law, at high energies the flattening 
of the spectrum due to Coulomb losses can be ignored, justifying the use of the ap- 
proximate single power-law form. The combined (primary plus secondary) electron 
spectral density distribution is then 

N,{y) = Ne.o{l+x)y-'^, (1) 
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Table 1 Star-forming galaxies: the data. 
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(Mpc) 


(kpc) 


(Jy) 


(cm^ 


3) (erg/s) (Ms/yr) (yr-') 


(Me,) 


(erg/s) 




Arp220 


74.7 


0.25 


0.3 


0.65 300 


45.75 


253 


3.5 


9 24+0"' 


< 42.25 


SB 


M82 


3.4 


0.23 


10.0 


0.71 200 


44.26 


8.2 


0.25 


9.37+!i;» 


40 21+""' 


SB 


NGC253 


2.5 


0.20 


5.6 


0.75 400 


44.23 


7.7 


0.12 


9.20+;i|Y 


39.76+«:|^ 


SB 


Milky Way 


- 


4.4 


- 


- 0.01 


43.75 


2.5 


0.02 


QO|+0.12 


38.91+«:i^ 


quiescent 


M31 


0.78 


4.5 


4.0 


0.88 0.01 


42.98 


0.43 


0.01 


9.88+0]^ 


38 66+"''' 


quiescent 


M33 


0.85 


2.79 


3.30 


0.95 0.03 


42.68 


0.22 


0.003 


Q lS+0'3 


< 38.54 


quiescent 


LMC 


0.049 


2.4 


285.0 0.84 0.01 


42.45 


0.13 


0.002 


8.86+0]2 


37.67+":^^ 


quiescent 


SMC 


0.061 


1.53 


45.3 


0.85 0.01 


41.45 


0.01 


0.001 


o fifi+0.03 


37-04+«:Ji 


quiescent 


NGC4945 


3.7 


0.22 


5.5 


0.57 300 


44.02 


4.7 


0.1-0.5 9.64+0^0 


40.30+«:i^ 


SB-FSy2 


NGC 1068 


16.7 


1.18 


6.6 


0.75 300 


45.05 


50 


0.2-0.4 9.71+[!]J 


41.32+«t^ 


SB4-Sy2 



I'lDistance (from Ackermann et al. 2012). 

PI Effective radius of star-forming region. See text. Data are from Persic & Rephaeli 2010 and 
refs. therein (Aip220, M82, NGC 253), Beck & Grave 1982 (M31), Tabatabaei et al. 2007 
(M 33), Weinberg & Nikolaev 2001 (LMC), Wilke et al. 2003 (SMC), Moorwood & Oliva 1994 
(NGC 4945), Spinoglio et al. 2005 (NGC 1068). 

1^1 1 GHz flux density Data are from Persic & Rephaeli 2010 and refs. therein (Arp 220, M 82, 
NGC 253)), Beck & Grave 1982 (M 31), Tabatabaei et al. 2007 (M 33), Klein et al. 1989 (LMC), 
Haynes et al. 1991 (SMC), Elmouttie et al. 1997 (NGC 4945), Kuhr et al. 1981 (NGC 1068). 

"iNon-thermal spectral radio index. Data are from Persic & Rephaeli 2010 and refs. therein 
(Arp 220, M 82, NGC 253), Beck & Grave 1982 (M 3 1), Tabatabaei et al. 2007 (M 33), Klein et 
al. 1989 (LMC), Haynes et al. 1991 (SMC), Elmouttie et al. 1997 (NGC 4945), Kuhr et al. 1981 
(NGC 1068). 

'^'Thermal electron density. Data are from Roy et al. 2010 (Arp 220), Petuchowski et al. 1994 
(M 82),Kewley et al. 2000 and Corral et al. 1994 (NGC 253), Cox 2005 (Milky Way), Beck 2000 
(M 31), Tabatabaei et al. 2008 (M 33), Points et al. 2001 (LMC), Sasaki et al. 2002 (SMC), Spoon 
et al. 2000 (NGC 4945), Kewley et al. 2000 (NGC 1068). 

^Total IR [i.e., (8 - 1000)^m] luminosity, in log (from Ackermann et al. 2012). 

PlStar formation rate, from SFR = Lir/(2.2 x lO''^ erg/s) (Kennicutt 1998). 

1*1 Core-collapse SN rate. Data are from Persic & Rephaeli 2010 and references therein (Arp 220, 
M 82, NGC 253), Diehl et al. 2006 (Milky Way), van den Bergh & Tammann 1991 (M 31, M 33, 
SMC, LMC; see also Pavlidou & Fields 2001), Lenain et al. 2010 and references therein 
(NGC 4945, NGC 1068). For NGC 1068 we also computed an upper limit to the SN rate 
(VsN < 0.39) using Mannucci et al.'s (2003) formula 

VsN = (2.4±0.1)xl0-2[LFiR/(10'"L,.)]yr-', being /fir = 1.26 X 10-11(2.58/60 -h/ioo) 
ergcm-^s-i (see Helou et al. 1988) with /go ^ 190Jy and /loo - 277 Jy. 

PI Gas mass (neutral plus molecular hydrogen: Mhi +Mh,), in log. Data are from: Torres 2004 for 
Arp 220; Abdo et al. 2010a for M 82, NGC 253, and the Milky Way; Abdo et al. 2010b for M 3 1 
and M 33; Abdo et al. 2010c for the LMC; Abdo et al. 2010d for the SMC; and Lenain et al. 2010 
for NGC 4945 and NGC 1068. 

''"iHigh-energy (>100MeV) y-ray luminosity, in log (from Ackermann et al. 2012). 
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where the electron Lorentz factor yis in the range 7i < 7 < 72, A^e.o is a normalization 
factor of the primary electrons, % is the secondary-to-primary electron ratio, and 
0' > 2 is the spectral index. Ignoring the contribution of low-energy electrons with 
7< 7i, the electron energy density is U^ = A^e.o (1 -^X)ni,tC^ [I- 7'^''d7, where 72 is 
an upper cutoff whose exact value is irrelevant in the limit of interest, 72 > > 7i . For 
q>2 and 72 >> 71, 



U,-Ne.Q{l+X)mecV, V{q-2). 



(2) 



For a population of electrons (specified by Eq. [U traversing a homogeneous mag- 
netic field of strength B in a region with (a spherically equivalent) radius rs located 
at a distance d, standard synchrotron relations yield 



A^e,o(l+Z) = 5.12^10-'' \ifa{q)-'B 



250 



9-1 



(3) 



where the scaled flux is /ighz Jy, a{q) is defined and tabulated in, e.g.. Tucker 
(1975), and v/= (i^)-^(i4)2(^). Use of Eq.© then yields 



f/e 



2.96 



X 10^2^2502 1^ 



ri 



-q+2 



{q-2)a{q) 



B 



ji+i 



(4) 



In order to compute f/e from Eq.© we need to specify 71 and B. To do so we 
make the following assumptions: 



(i) The low-energy limit of the electron power-law spectrum, 71 , marks the transition 
(for decreasing energy) from Coulomb (Rephaeli 1979) to synchrotron losses. For 
an electron of energy 7, the synchrotron loss rate is 



(^Y\ 



\ (Jf / syn 

whereas the Coulomb loss rate is 



= -1.30X 10 



-2\yif 



B 

\iiG 



„-i 






1.2 X \0-^\ 



e.th 



1.0- 



In (7/«e,th) 



75 



,-1 



(5) 



(6) 



(Rephaeli 1979). We then simply assume that electrons lose energy via Coulomb 
scattering for 7 < 7i and via synchrotron cooling for 7 > 7i . 

(ii) The particle energy density is in equipartition with that of the magnetic field, 
t^p + f^e = B^ l%%. In terms of the p/e energy density ratio, K, the equipartition con- 
dition is C/p[l + (1 +;t)/K-] = B'^jZn, so that 



B = 



7.44 X 10 
1+Z 



-21 



[1 



K- , ^"''250'?/V 



1+;^^ [q-2)a{q) 



5+9 



(7) 
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Inserting Eq.© into Eq.© we get (■gj-)syn °= 7i ' '' ■ Once the value of «e,th is 
specified (see Table 1), by equating Eqs.(|5]|6l) we deduce yi- 

The secondary-to-primary electron ratio j, which appears in Eq.©, depends on 
the injection p/e number ratio, rp/g = (mp/me)''''"J ' , and on the gas optical thick- 
ness to p-p interactions. Given the branching ratios in p-p collisions, only a third of 
these collisions produce electrons. The mean free path of CR protons due to p-p 
interactions in a medium of density n^ is App = (cTpp/ip)^'; for protons with kinetic 
energy T ~ few TeV the cross section is app ~ 50mb = 5 x 10^'^ cm^ (Baltrusaitis 
et al. 1984). For a typical SB ambient gas density «p ~ 150cm^^, App ^ 43kpc. 
The probability for a single CR proton to undergo a pp interaction in its 3D ran- 
dom walk through a region of radius r^ ^ 0.25 kpc (also typical of SB nuclei) is 
then V3'"s/App ~ 0.01. Thus, in a typical SB environment, characterized by rela- 
tively strong non-relativistic shocks (^inj — 2.2), the secondary to primary electron 
ratio is X = ;i^o V3 ('"s/App) ~ 0.3. In a more quiescent environment, with typical 
values /ip ~ 1 cm^^ and r^, ^ 2.5 kpc, x — 0.03. The higher value found in SBs is 
in approximate agreement with results of detailed numerical starburst models for 
energies ^ lOMeV (plotted in, e.g., Paglione et al. 1996, Torres 2004, De Cea et al. 
2009, Rephaeli et al. 2010). 

To compute the p/e energy density, K, we assume power-law spectra: (i) The 
electron spectral index q^ is deduced from the measured radio index a, generally 
qe = 2a + l. (ii) The proton spectral index is assumed to be close to the injection 
value, q^ ^ qiaj ~ 2.1 — 2.2, for the dense SB environments hosted in the central 
regions of some galaxies, and equal to the leaky-box value, q^ ~ qi„j + 5 ~ 2.7 
(where 5 ~ 0.5 is the diffusion index) for more quietly star forming galaxies. 

Finally, we obtain an explicit expression for Up: 



^^^i 



1+Z 



-1 



7.44x10-21 K , ^ ^250'5'/V 



[1 



1+Z ' 1+Z' iq~2)a{q) 



4 



^1 



(8) 



Using Eq.®, values of Up can be obtained from the relevant observational quan- 
tities for our sample galaxies; these values are listed in Table 2. The quantities in 
Eq.® are usually well determined for our sample galaxies, except for the p/e en- 
ergy density ratio k, for which a spectral index, q^, must be assumed. Given its 
possible values (i.e., (^p ~ 2.1 — 2.2 in SB regions, and (^p ~ 2.1 — 2.2 for quiescent 
galaxies), the uncertainty in the spectral index, 5qp ~ 0.1, translates to a factor of 
^2 uncertainty on K, i.e. typically an uncertainty of ^50% on Up as deduced from 
Eq.®. 



3.2 Energetic particles and y-ray emission 

Based on the calculation of 7-ray emission from SFGs outlined in Section 1, Up can 
be estimated directly from recent measurements of the nearby galaxies. In SBGs, 
such as M82 and NGC253, the central SB region (referred to also as the source 
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region) with a radius of ^ 200 — 300 pc is identified as the main site of particle 
acceleration. The injected particle spectrum is assumed to have an index q = 2, 
the theoretically predicted Np/Ns ratio is adopted, and equipartition is assumed. A 
measured radio index a ~ 0.7 in the source region implies ^ = 2a + l~2.4 there, 
indicating a substantial steepening due to diffusion {D °^ J )^ that cause the steady- 
state particle spectral index to he qo + 5 above some break energy. The procedure is 
similar when star formation does not (largely) occur in a burst in the nuclear region, 
but proceeds more more uniformly across the disk. 

For a source with ambient gas number density «gas, proton energy density Up, 
and volume V, the integrated hadronic emission from pp-induced Tt" decay is 

L>E = J 8>E'h-sUpdV s-\ (9) 

with the integral emissivity g^'^ in units of photon s^ ' (H-atom)^ '(e V/cm^)^ ' (Drury 
et al. 1994). Thus, Up can be determined from measurements of L>£ and ngas('"), and 
the particles steady-state energy distributions can be numerically calculated in the 
context of the convection-diffusion model. 

In addition to the high-energy detections of the two local SBGs M 82 and 
NGC253 (Abdo et al. 2010a, Acciari et al. 2009; Acero et al. 2009), several galax- 
ies with low SFR were also detected by the Fermi telescope, (i) the Andromeda 
galaxy M31 (Abdo et al. 2010b), with Up ~ 0.35 eVcm"^; (ii) the Large Mag- 
ellanic Cloud (LMC) whose average spectrum, either including or excluding the 
bright star-forming region of 30 Doradus, suggests Up ~ 0.2 — 0.3 eV cm^^ (Abdo 
et al. 2010c); (Hi) SMC for which Up ~ O.lSeVcm"^ was deduced (Abdo et al. 
2010d). For the Milky Way, the modeling of the Galactic diffuse HE emission along 
the lines outlined above requires an average C/p ~ 1 eVcm^^ (Strong et al. 2010; 
Ackermann et al. 201 1). Our values for Up determined from the measured GeV-TeV 
fluxes are listed in Table 2. 



3.3 Energetic particles and Supernovae 

The SN origin of energetic particles suggested early on; as a test of this hypothesis, 
we estimate of Up by combining the SN rate with the proton residence time, Tres. 
assuming a fiducial value for the fraction of SN kinetic energy that is channeled 
to particle acceleration. The residence time is determined from the p-p interaction 
time, and the two propagation timescales of advection and diffusion; 
f;j The energy-loss timescale for/?/? interactions, Tpp = (appCWp)"'; for protons with 
kinetic energy E^ 10 TeV for which C7pp ~ 50 mb, this timescale is 

Tpp - 2xl0^( "P )-'yr. (10) 

MOOcm ^' 
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(ii) Particle advection out of the disk mid-plane region in a fast SB-driven wind oc- 
curs on a timescale Tadv determined from the advection velocity for which we adopt 
(except where noted otherwise) the nominal value Vadv ^ 1000 km s^^ deduced 
from measurements of the terminal outflow velocity of ~^ 1600-2200 km s^' in M 82 
(Strickland & Heckman 2009; see also Chevalier & Clegg 1985 and Seaquist & Ode- 
gard 1991). For a homogeneous distribution of SNe within the SB nucleus of radius 
Ts, the advection timescale is 

Tadv=i7.5xl0n— ^^) ( ^21!^)-iyr. (11) 

^OJkpc^ MOOOkms-i^ ^ 

(Hi) As noted in the previous section, diffusion is likely to be random walk against 
magnetic field inhomogeneities, with an estimated central diffusion coefficient, D ~ 
3 X 10^*^ cm^/s, assuming a magnetic coherence scale of A ~ 1 pc. Thus, diffusion 
out of the central 0.5 kpc is estimated to occur on a timescale 

Now, since the weighted residence time is 

"^res = "^pp + "^adv + ''■diff ' (^^) 

it is expected that - under typical conditions in central SB regions - p-p collisions 
and advection, more so than diffusion, effectively determine the survival there of 
energetic protons. 

During Tres, the number of SN is VsN^res; the kinetic energy deposited by each of 
these into the ISM is Egj = 10^' erg (Woosley & Weaver 1995). Arguments based 
on the energetic particle energy budget in the Galaxy and SN statistics suggest that 
a fraction 77 ^ 0.05 — 0.1 of this energy is available for accelerating particles (e.g., 
Higdon et al. 1998). Thus, the proton energy density can be expressed as 

The resulting values of Up in the sample galaxies are listed in Table 2. 



4 Discussion 

Detections of VHE 7-ray emission associated with ongoing star formation in M 82 
and NGC 253 add significant new insight on the enhanced energetic electron and 
proton contents in SBGs, and on their propagation in disks of spiral galaxies. The 
common star-driven nature of NT phenomena in the wide class of non-AGN SFGs 
implies that energetic particle densities and radiation fields are self-similarly scaled 
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Table 2 Star-forming galaxies: Proton energy densities+. 



Object 


7-ray 


radio 


SN 


other 


'•s T,-es 




meth. 


meth. 


meth. 


meth. 


(kpc) (yr) 


Arp220 


- 


1027 


1142 


_ 


0.25 2.0E+4 


M82 


250"''= 


250 


234 


- 


0.23 4.5E+4 


NGC253 


220''''= 


230 


213 


- 


0.20 6.7E+4 


Milky Way 


Id 


- 


1.2 


IJ 


4.4 2.0E+7 




6'= 


- 


5 


- 


0.2 2.5E+6 


M31 


0.36f 


0.22 


0.7 


- 


4.5 2.5E+7 


M33 


< 0.43^ 


0.38 


0.7 


- 


2.8 2.0E+7 


LMC 


0.21-0.318 


0.22 


0.4 


- 


2.5 4.4E+7 


SMC 


O.IS'' 


0.39 


1.1 


- 


1.5 1.4E+7 


NGC4945 


200' 


201 


215 


- 


0.22 4.5E+4 


NGC 1068 


- 


65 


61 


- 


1.2 l.OE+6 


+ Values ai-e in eVcm^^. 











(a) Acciari et al. (2009; see also Persic et al. 2008, De Cea et al. 2009). (b) Acero et al. (2009). (c) 
Abdo et al. (2010a). (d) Strong et al. (2010). (e) Aharonian et al. (2006). (f) Abdo et al. (2010b), 
and Drury et al. (1994) for M 33. (g) Abdo et al. (2010c). (h) Abdo et al. (2010d). (() Lenain et al. 
(2010). (j) Webber (1987). 



with SF activity, from quiescent systems to intense SBGs, in spite of the wide range 
of intrinsic physical conditions in these systems. 

We have briefly outUned our treatment of the steady state spectro-spatial distribu- 
tions of energetic electrons and protons in SFGs, exemplified in the case of the two 
nearby SBGs M82 and NGC 253. This approach is based on a numerical solution 
of the diffusion-convection equation for particle distribution functions, following 
the evolution from the acceleration sites throughout the disk as the particles lose 
energy and propagate outward. Key observational normalization is based on mea- 
surements of radio synchrotron emission which, through an initial (theoretically 
assumed) Np/Ng ratio provides also the normalization of the proton component. As- 
suming equipartition then allows to relate the local value of the mean magnetic field 
to the particle energy densities. The numerical solution of the diffusion-convection 
equation is based on an iterative procedure to determine the particle densities and 
mean field strength in the galactic center, and evolve these quantities by accounting 
for all relevant energy losses, and normalizing central values of these quantities by 
fitting to the measured radio spectrum from the central galactic (or SB) region. The 
quantitative viability of this approach is confirmed by the good agreement between 
the predicted high energy emission from M 82 and NGC 253 and measuremnts with 
Fermi, H.E.S.S., and VERITAS. 

Significant detections of the NT emission from above two SBGs at lower (below 
100 keV) X-ray energies, as would be expected by the currently operational NuS- 
TAR telescope, will provide additional spectral coverage that will allow separating 
out the spectral electron and proton components. NuSTAR is the first X-ray tele- 
scope with capability to resolve this emission; if this is indeed achieved, important 
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new insight will be gained on the evolution of the electron spectro-spatial distribu- 
tion across the disks of these nearby SBGs. 

The three methods we have discussed to estimate energetic particle energy densi- 
ties are clearly not independent. The y-ray method and the radio method are coupled 
through the p/e ratio at injection, through the secondary-to-primary electron ratio, 
and through the imposed condition of particle-field equipartition. The SN method 
is not independent of the 7-ray method either, because both depend on the proton 
residence time, although - unlike the 7-ray and radio methods - it does not depend 
on the particle radiative yields but on the statistics of core-collapse SN. Also, the 
three methods do not stand on equal footing: with the 7-ray, radio, and SN methods 
we, respectively, either measure, infer, and or estimate the value of Up. A substan- 
tial agreement among estimates based on the three methods is found for most of 
the galaxies in Table 1. The only exceptions are the SMC and NGC 1068. As for 
the former, the proton confinement volume could be small, so that most particles 
diffuse out to intergalactic space (Abdo et al. 2010d). If so, the 7-ray method yields 
the (lower) actual proton energy density, whereas the radio and SN methods esti- 
mate the (higher) produced amount. NGC 1068 hosts a prototypical Seyfert-2 nu- 
cleus (e.g., Wilson & Ulvestad 1982) surrounded by a spherical circumnuclear SB 
shell with external radius of 1.5 kpc and thickness 0.3 kpc, and mass 3.4 x IO^Mq 
(Spinoglio et al. 2005); its implied energy density is Up « 65eVcm^^ from both 
the radio and SN methods. [We note that Lenain et al. (2010) suggested that the HE 
emissions of the above two SB and Seyll galaxies NGC 4945 and NGC 1068 are 
powered by, respectively, star formation and AGN activity.] 

A debated aspect of proton energy loss and propagation times is whether the for- 
mer is shorter than the latter; if so, the system is said to be a 'proton calorimeter' 
(e.g., Lacki et al. 2010, 2011). No galaxy in the above sample is found to be in 
the calorimetric limit; the two SBGs M82 and NGC 253 would seem to be only 
marginally close to this limit. The presence of fast, SB-driven galactic winds ad- 
vecting energetic particles out of the disk seems to be a ubiquitous feature in SBGs, 
limiting the degree at which they can be calorimetric. More generally, it is known 
that energetic particles do diffuse out of non-AGN SFGs, as evidenced also by sig- 
nificant radio emission (Ferrari et al. 2008), and possibly also high energy NT X-ray 
emission (Rephaeli et al. 2008) from large central regions of galaxy clusters. The 
estimated intracluster particle (and indeed also the magnetic) energy densities are 
sufficiently high, suggesting origin in the cluster galaxies. 

Based on the reasonable hypothesis that local SB galaxies resemble young galax- 
ies which were particularly abundant in the early universe, their contributions to the 
X-7-ray backgrounds are of obvious interest (e.g., Rephaeli et al. 1991). Calcula- 
tions of the superposed emission from SBGs (Pavlidou & Fields 2002, Persic & 
Rephaeli 2003, Thompson et al. 2007, Fields et al. 2010, Lacki et al. 2011, Steck- 
ers & Venters 201 1) indicate that this emission constitutes at least a tenth of these 
backgrounds. 
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